Co-tunneling current and shot noise in quantum dots 
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We derive general expressions for the current and shot noise, taking into account non-Markovian 
memory effects. In generalization of previous approaches our theory is valid for arbitrary Coulomb 
interaction and coupling strength and is applicable to quantum dots and more complex systems 
like molecules. A diagrammatic expansion up to second-order in the coupling strength, taking 
into account co-tunncling processes, allows for a study of transport in a regime relevant to many 
experiments. As an example, we consider a single-level quantum dot, focusing on the Coulomb- 
blockade regime. We find super-Poissonian shot noise due to spin-flip co-tunneling processes at an 
energy scale different from the one expected from first-order calculations, with a sensitive dependence 
on the coupling strength. 

PACS numbers: 73.63.-b, 73.23.Hk, 72.70.+m 



Introduction. - The study of shot noise in trans- 
port through mesoscopic devices, such as quantum dots 
or molecules, has become a field of intense theoretical 
and experimental research^. It provides additional in- 
formation, not contained in the current, about system 
parameters that govern the electronic transport |2( • For 
weak coupling between the dot and metallic electrodes, 
transport is dominated by sequential tunneling, usually 
described by first-order perturbation theory in the cou- 
pling strength. Shot noise in this limit has been stud- 
ied for a variety of models, including effects from elec- 
tronic interactional^, multi-level dots0,|f|, inclusion of 
photonic and vibrational modes @,@], or spin-polarized 
leads0- However, a theory taking fully into account in- 
teraction effects as well as arbitrary coupling strength 
has not been presented. Expressions for the shot noise 
in terms of non-equilibrium Green functions have been 
derived only in the absence of Coulomb interactions 8] or 
in a perturbative expansion thereof @, [13 ■ 

Shot noise has been measured in various experimen- 
tal realizations [HI Il3| . As the strength of the cou- 
pling in a given experiment is not a priori known, it 
is unclear whether first order calculations are sufficient. 
Second-order tunneling (co-tunneling^ij]) processes can 
play an important role for the conductance particu- 
larly in the Coulomb-blockade regime. Shot noise in this 
regime has been discussed in Ref.llal where the possibility 
of super-Poissonian noise was suggested. However, that 
work was limited to pure co-tunneling processes, which 
is too restrictive for the experiment of Ref. flil 

In this Letter we provide a general formulation of 
the current and shot noise, valid for arbitrary coupling 
strength, while accounting fully for the Coulomb interac- 
tion. By making use of a diagrammatic formulation we 
expand the expressions up to second order in the coupling 
strength and study as an example a single-level quan- 
tum dot, which is out of equilibrium due to an applied 



bias voltage. The theory covers the Coulomb-blockade 
regime (low bias), the sequential tunneling regime (bias 
larger than the first single-charge excitation energy) and 
the crossover between both. We find super-Poissonian 
shot noise in the Coulomb-blockade regime due to spin- 
flip co-tunneling processes at a different energy scale than 
expected from first-order calculations. Furthermore, we 
show that the noise to current ratio is highly sensitive to 
the coupling strength. This may serve as an additional 
spectroscopic tool for the couplings. 

The model. - The Anderson impurity model is 
based on the Hamiltonian H = Hi + tin + ti^ot + 
Ht,l + Ht,r with ti r = J2ka e ka r a\. ar ak a r, H^ot = 



4, 



Un-\ni and H^ r = J2ka(^r a 



l kcrr l 



h.c.) 



for r = L,R. Here, Hi and Hn describe the left and 
right electrode with non-interacting electrons, -ffdot the 
quantum dot with one (spin-dependent) level of energy 
e a and Coulomb interaction U for double occupancy of 
the dot. Tunneling between the leads and the dot is mod- 
eled by Hi i and -£/t,r- The coupling strength is char- 
acterized by the intrinsic line width T r = 2itp e \t r \ 2 , and 
r = Ti + Tn, where p e is the (constant) density of states 
of the leads. The creation operators aL and c£ refer to 
the electrodes and the impurity (dot), and n a = c\c a . 

We are interested in the current I and the (zero- 
frequency) current noise S, which for eV 3> k^T is 
dominated by shot noise. Both are related to the (sym- 
metrized) current operator I = (In — Ii)/2 with I r = 

from which we obtain 



*(eA)E feCT ( 
I = (I) and 
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dt 
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Diagrammatic technique. - In Ref.0we formulated 
a theory of current noise for transport to first order in 
r, based on a diagrammatic language that has been de- 
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velopped for a systematic perturbation expansion of the 
current through localized levelsfl?). All transport prop- 
erties are governed by the non-equilibrium time evolution 
of the density matrix. The electrode degrees of freedom 
can be integrated out, and we obtain a reduced density 
matrix for the dot degrees of freedom only, labeled by 
X- The time evolution of the reduced density matrix, de- 
scribed by the propagator H x i x (t', t) for the propagation 
from a state x at ti me t to a state x' at time t', can be vi- 
sualized by diagrams 0,^3] on the Keldysh contour. The 
full propagation is expressed as a sequence of irreducible 
blocks (self energies) W x > x (t',t) that are associated with 
transitions from state x a t time t to state x' & t time t' . 
This yields the Dyson equation 



/t rti 
dt 2 J diiW(t 2 ,*i)II(ti,t) 



(2) 



for the propagator, where the bold face indicates ma- 
trix notation related to the dot state labels. For the 
following, it is convenient to introduce the Laplace trans- 
form W(z) = ft f dt e zt W(0, t) and the definitions 
W = W(z)\ z=0+ and°<9W = (dW(z)/dz) | z=0+ . 

In the long-time limit, i.e. for time differences t' — t 
larger than the correlation time (oc 1/T) over which the 
system forgets its initial state, the propagator becomes 
p st ® e T , where e T = (1, . . . , 1), and p st is the vector of 
the stationary probabilities, determined from Wp st = 0, 
independent of the initial (diagonal) density matrix. As 
W has a zero eigenvalue, it cannot be inverted. With 
the normalization condition e T p st = 1 we obtain the 
stationary probabilities p st by solving 



(w P st j 



(3) 



where W is identical to W but with one (arbitrarily cho- 
sen) row xo being replaced by (r, . . . ,r) L 2j. For a dia- 
grammatic representation of the current, we introduce a 
block W 7 , in which one (internal) vertex is replaced by 
an external due to the current operator. We get 



2ft v 



(4) 



The shot noise, Eq. (JTJ, involves expectation values 
of two current operators. They either appear both in a 
single irreducible block, which we denote by W 77 , or in 
two different blocks W 7 . We find 



e T W = 0. The set of matrix equations Eqs. J2J - © 
constitute the starting point for all numerical results pre- 
sented below. For a systematic perturbation expansion 
of current / and noise S in T, we expand all quantities 
W, W 7 , W 77 , aW, aW 7 , p st and P order by order. 
We remark that for transport in first order (sequential 
tunneling), the above expressions simplify, as all contri- 
butions involving the derivatives 5W and <9W 7 disap- 
pear as a consequence of the fact that W starts at order 
r, p st at r , and P at T -1 . These derivatives are asso- 
ciated with non-Markovian behavior of the system, and 
have not been taken into account in Refs. ©andi They 
are absent for first-order transport but are important for 
second- and higher-order corrections. (For a discussion 
of non-Markovian effects see also Ref.[lJJ) Higher deriva- 
tives will not appear for the shot noise even for higher- 
order corrections. 

In the following we will focus on the Anderson impurity 
model, but we emphasize that the expressions presented 
above also cover more general situations, such as quan- 
tum dots or molecules with many levels and more general 
forms of the interactionsll9, 201. 
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FIG. 1: Current I and shot noise 5* vs bias voltage for feT = 
O.lmeV, e L = -1.5meV, e T = 0.5meV, U = 0, and T L = Tr = 
r/2. First order (dotted lines) and second order (dashed lines) 
are compared to the exact results (solid lines) for F /k&T = 



S = —e T [W 77 + W 7 (PW 7 + p st ® e T <9W 7 )] p st 

(5) 

with <9W 7 = (dW I (z)/dz) \ z=0+ . The object P = 
J°oo dt j- [11(0, t) - n(0, -oo)] is determined by 

WP = (p st $ e T - 1)(1 - 5 X ,, X0 ) - 5Wp st ® e T . (6) 

Here we use the extra condition e T P = 0, which fol- 
lows from the definition of P, the Dyson equation, and 



Results. - We recall first the main features of 
first-order (sequential) transport. Both current and 
noise increase monotonically with bias voltage displaying 
plateaus separated by thermally broadened steps. The 
step positions are determined by the energy of single- 
charge excitations. The coupling parameters set the 
plateau heights. 

Higher-order processes modify the current and noise 
in two different ways. First, they increase the width of 



3 



the steps, which effectively is given by the sum of T and 
T. Second, they allow for transport in the Coulomb- 
blockade region at low bias, where sequential tunneling is 
suppressed. With increasing coupling strength T, second- 
and eventually also higher-order corrections to transport 
become more and more important. To illustrate the va- 
lidity range of our second-order perturbation expansion 
we first consider the non-interacting limit, U = 0, since 
in this case exact results 0, are available for the cur- 
rent I u=0 = e/h J ^E/jHI/l^) -/rM) an d shot 
noise S u=0 = 2e 2 /hJ du^M^lfL^O- ~ + 

/HM(1 " /rM)] + T CT H(1 - T a (L0))[f L (L0) - / R (u,)] 2 } 

with Ta (u) = r L r R /[(Lo - e a ) 2 + (T/2) 2 ]. In Fig. □ we 
compare current and noise in first and second-order with 
the exact results. We choose a symmetric bias voltage 
(/jl = —Mr = eV/2). Outside the Coulomb-blockade 
regime, the second-order corrections (dashed lines) to 
sequential tunneling (dotted lines) start to become im- 
portant for T/k B T < 0.5. As long as T/k B T < 1, the 
exact curves (solid lines) are perfectly reproduced by 
second-order perturbation theory, while the sequential- 
tunneling results clearly deviate. For T/k B T = 2 third- 
order contributions start to play a role, at least for the 
noise, where unphysical non-monotonicities arise around 
the steps. We, therefore, restrict ourselves in the follow- 
ing discussion to T < k B T. 
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FIG. 2: Current I and conductance dl/dV (inset) vs bias- 
voltage for ej = — 1.5meV,e| = 0.5meV, U = 4meV and 
Tl = Tr = T/2 for various values of F and ksT. The broad- 
ening of the first step due to F and ksT is shown in the inset. 
The dashed and dotted curves with the same (F + /cbT) have 
about the same width. 

The elastic co-tunneling processes which do not change 
the dot state or its energy allow for an electron exchange 
with the reservoirs via an intermediate virtual state. This 
leads to a finite linear conductance G — dl /dV\v=a- 
Also the noise persists at zero bias, S — 4k B TG, consis- 
tent with the equilibrium fluctuation-dissipation theorem 
(FDT). In the Coulomb blockade regime the FDT can 



be extended to non-equilibrium 0] and takes the form 
S {2) {V)/2eI {2) {V) = coth(eV/2k B T). Our theory satis- 
fies this relation, however, we stress that it holds only 
in the regime of purely elastic co-tunneling processes. In 
Fig. |21 we show the current / normalized to ireT/h for 
the same set of energy parameters as in Fig. ^but with 
a finite interaction U = 4meV. Since the bias is ap- 
plied symmetrically, the dot preferably occupies the state 
with spin J, (Coulomb-blockade) until it can be emptied 
due to first order hopping processes around 3mV (first 
step). Further steps arise around 5mV and 9mV due 
to the double occupied dot state. This parameter set 
is similar to the experimental situation of Ref. for a 
quantum dot with occupation N = 2 21]. In Fig. 3 of 
that paper, a conductance step is observed inside the 
Coulomb blockade diamond, that is attributed to inelas- 
tic co-tunneling processes. For our model one expects 
this inelastic co-tunneling feature in the conductance at 
a bias of e co /e — (ef — ej.)/e = 2mV. This feature is 
hardly noticable in the conductance plot of the inset in 
Fig. |21 because the coupling V is relatively weak and the 
energy e co is fairly close to the sequential-tunneling en- 
ergy. However, the inelastic co-tunneling processes can 
clearly be observed in the shot noise and the Fano factor 
F = S/2eI discussed below. 

We note that the dashed and dotted curves in Fig. [21 
with same total sum (r + k B T) almost lie on top of 
each other. The differential conductance plot (inset) 
shows that the temperature effect is a little stronger: the 
dashed curve with the highest temperature has the lowest 
peak. The FWHM of the conductance peaks is between 
0.5mV - 0.8mV - 6(r + /cB^has compared to 5.Uk B T 
for pure sequential tunneling 15]. We also note a shift 
of the peak position from the sequential-tunneling value 
3mV to somewhat lower bias voltages, indicating a renor- 
malization of the level positions. 

The Fano factor F = S/2eI for a fixed tempera- 
ture and a sequence of coupling constants T, covering 
three orders of magnitudes, is shown in Fig. At low 
bias, the Fano factor varies [13 as coth( e y/2A: j BT) un- 
til it reaches the value 1, as expected for uncorrelated 
systems. For bias voltages around the spin-flip excita- 
tion energy e rn = 2meV, the Fano factor becomes super- 
Poissonian 16], F > 1. Once sequential tunneling be- 
comes dominant (at a bias > 3mV), it drops to values 
between 1 and 1/2. 

The super-Poissonian Fano factor appears for bias volt- 
ages at which the spin-j level acquires some finite oc- 
cupation probability, either due to inelastic spin-flip co- 
tunneling, or due to sequential tunneling processes that 
are exponentially suppressed but, for the chosen param- 
eters, still finite. The enhancement of the noise is mainly 
described by the second and third term of Eq. 10). It 
describes bunching of the transfered t-electrons during 
the time when this transport channel is not blocked by 
the dot being occupied with a J,-electron. Spin flip pro- 
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FIG. 3: Fano factor F — S/2eI vs bias voltage for the same 
parameter set as Fig.[2]but fixed temperature k&T — O.lmeV 
and various V. Inelastic co-tunneling leads to a super- 
Poissonian Fano factor if the bias is larger than e co /e = 2mV. 
First order processes may also lead to a super-Poissonian 
value at a scale 2tj /e = lmV. The crossover between these 
energy scales runs over three orders of magnitude in the cou- 
pling. Outside the Coulomb blockade regime the first order 
results are recovered already at Y/ksT ~ 0.1. The inset shows 
a sketch of the transport situation at eVbias = ^co- 
cesses leading to super-Poissonian noise have also been 
discussed in Ref. in the case of ferromagnetic leads. 
Both the position and the height of the peak in the Fano 
factor depends on all system parameters. In Fig. |31 we 
study the dependence on the ratio Y/k^T. 

With reduced coupling strength Y, the peak increases 
and moves towards lower bias. For Y/k^T = 0.0005 
our result (dotted line) coincides with that of a pure 
first-order calculation. At larger coupling, though, the 
second-order terms lead to a severe modification of the 
result because the rates in the Coulomb-blockade regime 
depend only algebraically and not exponentially on the 
energy. We emphasize that, since the peak is close to the 
onset of sequential tunneling, an analysis purely based on 
co-tunneling processes^ would not be sufficient either. 
The range of Y /k^T ratio over which second-order dom- 
inated co-tunneling crosses over to sequential tunneling 
is rather large, i.e., the Fano factor provides a sensitive 
measure of the coupling strength. For some range of the 
ratio Y /k-gT the inelastic co-tunneling effects on the shot 
noise can be measured experimentally [2^. The impor- 
tance of the second-order processes for the peak of the 
Fano factor contrasts with the situation at larger bias 
beyond the sequential tunneling threshold when second- 
order corrections only become noticable for Y /k^T ~ 0.1. 

In summary, we presented a theory of current and 
shot noise within a diagrammatic technique that includes 
higher-order tunneling processes in the coupling of a 
quantum dot to metallic electrodes. As an example, 
we studied current and noise for an Anderson impurity 
model with finite spin splitting. We showed that espe- 
cially the steps and the Coulomb-blockade regions are 



strongly affected by second-order processes and provide 
additional information compared to first order. Spin- 
flip processes of different origin in first and second order 
lead to super-Poissonian noise in the Coulomb-blockade 
regime. This provides a spectroscopic tool to characterize 
the dot-electrode couplings by a combined measurement 
of current and shot noise. 
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